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For a massless gas with constant cross section in a homogeneous, isotropically expanding spacetime 
we reformulate the relativistic Boltzmann equation as a set of non-linear coupled moment equations. 
For a particular initial condition this set can be solved exactly, yielding the first analytical solution of 
the Boltzmann equation for an expanding system. The non-equilibrium behavior of this relativistic 
gas can be mapped onto that of a homogeneous, static non-relativistic gas of Maxwell molecules. 

PACS numbers: 25.75-q,51.10.+y,52.27.Ny,98.80.-k 


1. Introduction. The relativistic Boltzmann equation 
plays a prominent role in understanding the complex 
non-equilibrium dynamics displayed by dilute relativis¬ 
tic gases. It has applications in many areas of physics 
including, e.g., theoretical description of the quark-gluon 
plasma [ll-@ , neutrino transport in supernovae 0,0 , and 
structure formation in cosmology |9l4llj. While analyti¬ 
cal solutions of the Boltzmann equation have been thor¬ 
oughly studied for homogeneous systems 12-[l4|, solu¬ 
tions for an expanding system remain to be found even 
in the non-relativistic regime. 

In the relativistic regime, progress in this direction was 
made recently |15l422l | using the Anderson-Witting equa¬ 
tion [ 2 ^ . an approximation of the Boltzmann equation 
that relies on the relaxation time approximation |24l | . In 
this scheme, the nonlinear collision kernel of the Boltz¬ 
mann equation is replaced by a linearized version that 
qualitatively describes the relaxation of the system to 
equilibrium on a microscopic time scale. These analyses 
served to improve our understanding of the domain of ap¬ 
plicability of a number of extended hydrodynamic theo¬ 
ries used in the description of rapidly expanding plasmas 
H-i^, with focus on the description of ultrarelativistic 
heavy-ion collisions. 


However, a complete description of dilute gases can 
only be achieved by solving the full Boltzmann equa¬ 
tion. While this can be done numerically, simple yet 
physically motivated analytical solutions of the Boltz¬ 
mann equation, if available, can lead to valuable insights 
into non-equilibrium phenomena. In this Letter we take 
a step in this direction and derive the first analytical so¬ 
lution of the full Boltzmann equation for an expanding 
dilute gas. This is done using the method of moments 
[1, [ 2 ^ to calculate the nonlinear collision term of the 
relativistic Boltzmann equation for a massless gas with 
constant cross section in a homogeneous and isotropically 
expanding spacetime. The derived solution describes the 
expansion-driven deviation of the particles’ momentum 


distribution from local thermal equilibrium as a func¬ 
tion of the system’s expansion rate. We note that the 
non-expanding limit of this solution, which describes the 
homogeneous relaxation of a relativistic gas towards equi¬ 
librium, is also new. 

2. Boltzmann equation. We consider a homogeneous and 
isotropically expanding system of massless particles in 
a Friedmann-Lemaitre-Robertson-Walker (FLRW) met¬ 
ric [ 2 ^ [13 (the most general homogeneous and isotropic 
metric in flat space) 

ds^ = dt^ — a^{t) {dx^ -I- dy^ -\- dz^) . (1) 

We note that in this paper we make the further restriction 
of zero global curvature. For a gas of massless particles, 
the energy density and the pressure are related as £ = 3p, 
and here the scale factor a(t)>0 is a free function 
that may be fixed by additional physics assumptions. 
For the FLRW metric the nonzero Christoffel symbols 
are Fq^ = <5* H{t) and = a{t)d{t) Sij, where i,j denote 
spatial indices, H{t) = d{t)/a{t) is the Hubble param¬ 
eter, and 1/—g = a^{t), with g being the determinant 
of the metric in ©• Even though the fluid flow of this 
system is locally static, = (1,0, 0,0), the expanding 
FLRW geometry induces a nonzero fluid expansion rate 
0{t) = = 3H{t). 

The dynamics of the single-particle distribution func¬ 
tion, fk{x), is given by the relativistic Boltzmann equa¬ 
tion in curved space [ll, [I^, [^sl - ls^l 

k>^d^h + r^^Mk>^^^=c[f]. ( 2 ) 

Since the FLRW geometry is based on the assumptions 
of spatial homogeneity and isotropy, fk cannot depend 
on spatial position x and must be locally isotropic in 
momentum, depending only on u-k = k^ where for mass¬ 
less particles = k/a{t) with k=\k\. We therefore write 
fk(x) = fk{t) from here on. 

The symmetries of the FLRW metric strongly con¬ 
strain the form of the energy-momentum tensor and 
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particle 4-current iV^ of the matter. Due to local momen¬ 
tum isotropy the viscous shear-stress tensor and particle 
diffusion current vanish exactly. Further, the bulk vis¬ 
cous pressure is zero for massless particles. Therefore, the 
conserved currents always take their equilibrium form, 
and the time evolution of the energy and particle den¬ 
sities, e and n, is fully determined by the conservation 
laws 


dtn + 3nH{t) = 0, dte + AeH^t) = 0. (3) 

With initial condition a{to) = l, they are solved by 
n{t) =n(to)/a^(t) and e{t) = e(tQ)/a^(t). 

Even though the conserved currents = eu^u'^ — 
pg^'^ and = nu^ take the same form as in local ther¬ 
mal equilibrium, the system itself does not have to be 
in equilibrium. In fact, the local momentum distribution 
fk{t) is driven away from its local equilibrium form by an 
amount proportional to the expansion rate 9{t) = 3H{t) 
of the FLRW geometry. However, this process happens 
without disturbing the spatial homogeneity and isotropy 
of the system and, consequently, there are no out-of- 
equilibrium corrections to the conserved currents. A sim¬ 
ilar behavior was found in [HI [HI in the context of the 
relaxation time approximation. 

The contributions from the Christoffel symbols in 
Eq. (12) cancel exactly, and the equation reduces to 

k°dtfk=C[f]. (4) 


Ignoring quantum statistics, the collision term reads 

C[f] = \ j Wuk'^pp' ifpfp'-fkfk'), (5) 

^ J k'pp' 

where = J [(27r)^-^—g and Wkk'^pp' is the 
transition rate. To make progress on evaluating m we 
make the simplifying assumption of isotropic scattering 
with energy-independent total cross section a. Then 
Wkk'^pp' takes the form [28l - IHj| 

Wkk'-^pp’ = sS^{k+k'-p-p), (6) 

where s = {k^+k'^){k^ + k'^). Thus, the Boltzmann equa¬ 
tion becomes 

k°dtfk = s6'^{k+k'-p-p'){fpfp,-fkfk')- 

^ J k'pp' 

( 7 ) 

Note that even for the highly symmetric case considered 
here the relativistic Boltzmann equation is still a non¬ 
linear integro-differential equation for fk ■ 

3. Moment equations. We proceed to solve the Boltz¬ 
mann equation using the method of moments 0. Due 
to local momentum isotropy, the distribution function fk 
can be fully described by scalar moments only [Hj |: 


Pm{t) = [ {u-kr+^ fk{t) = [ fkit) 

Jk Jk 


( 8 ) 


The time dependence of the two lowest moments, the 
number density n{t) = pQ{t) and the energy density 
e{t) = pift), is given by the conservation laws already 
discussed. For a classical gas of massless particles, they 
provide the time-dependence of the temperature T and 
fugacity A= exp(^/r) (where p. is the chemical poten¬ 
tial) via the matching conditions Q T = e/{3n) and 
Vd^ = jT'^ where Vd is the number of massless degrees 
of freedom (including statistical degeneracy factors). In 
the following we set = I for simplicity. 

From the matching conditions it follows that 
T{t)=T{to)/a{t) and p{t) = pfto)/a{t) such that the 
fugacity A is time-independent. T{t) and A define 
the local equilibrium distribution function /^'^(t) = 
Aexp(—u-fc/T(t)) and the equilibrium scalar moments 

= (9) 

To express the Boltzmann equation in terms of the scalar 
moments Pm we multiply © by {u-k)^ and integrate over 
k. This results in 

dtPn,{t) + {3+m)H{t)p^{t) = (10) 

where the gain and loss terms are defined by 

s {u-p)"^d'^{k+k'-p-p')fkfk', 

W = s {u-k)"^6^{k+k'-p-p')fkfk'- 

^ J kk'pp' 

( 11 ) 

Here the gain term was simplified by using the symmetry 
of the transition rate under interchange of the incoming 
and outgoing momenta. With 

f S'^ik+k'-p-p') = (12) 

Jpp> 

the loss term reduces straightforwardly to Q 

= ^Pm{t)n{t). (13) 

However, the gain term is more involved. One first writes 


cSl(.t) = ^[ sfkfk'V^, (14) 

where the inner kernel is 

Pm = (27r)^ / (u-p)"^ S'^ik+k'-p-p'). (15) 

Jpp' 

This quantity is a scalar and can be calculated in 
any frame. In the center of momentum frame, the 4- 
dimensional delta function can be integrated and the re¬ 
maining angular integrals can be performed analytically 
[ 3 ^ . The final result can be expressed as follows 


m! {k+k')'^+^-i \k+k'\>~^ 

[2a(0]’" {m+l-j)\ j\ 


(16) 


27r^ 


0 


dk k^~'~‘^ fk{t) (m e No). 
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Plugging this back into (fl3) and performing some alge¬ 
braic manipulations one obtains 


= ,17, 

Using these results in Eq. (flOl) we derive the following 
exact evolution equation for the moments: 


dtPm{t) + [{i+m)H{t) + an{t)]pm{t) 
Pj (t) Pm-j (t) 


= 2 (m-|-2) to! cr jr 


3=0 


{j+2)\ {m-j+2)1' 


(18) 


This infinite set of coupled nonlinear differential equa¬ 
tions for the scalar moments Pm{t) is equivalent to the 
original relativistic integro-differential Boltzmann equa¬ 
tion. For TO = 0 and to = 1 Eq. (flSl) reduces to Eqs. 
for the particle and energy densities. Higher-order mo¬ 
ments are more sensitive to the distribution function 
at higher momentum. We can further simplify the 
moment equations by introducing the scaled moments 
Mm{t) = pm{t)/p^(t) [ 3 ^ and the scaled time i = tjio 
where £0 = 1 /{a n{to)) is the (constant) mean free path 
at time to- This yields the surprisingly simple evolution 
equations 


\i) 


di 


+ Mm{i) — Mj{i)Mm-j{i)- (19) 


TO-l-1 


1=0 


Equation m is the main result of this Letter. While 
the nonlinear coupling between different moments was 
expected from the nonlinearity of the collision kernel, 
the same cannot be said about another key feature of 
Eq. (fT^ : it can be solved recursively, i.e. the solution 
of the evolution equation for requires only previ¬ 
ously solved moments Mk {t) of lower order k <n. This 
property depends on our choice of an energy-independent 
cross section; it is essential for being able to solve m 
analytically. 

All information about the expansion appears in the 
factor a^{t) multiplying the time derivative in (1191) . 
Whether or not local equilibrium can be achieved 
thus depends on the state of expansion of the sys¬ 
tem. For example, when a{i) the mean free path 

£{t) = l/{(Tn{t)) increases faster than the expansion rate, 
limj^j^ £(t)0(t) —>■ 00 , and local equilibrium cannot be 
reached even at asymptotically large times. Moreover, 
if the initial /fe(to) is positive definite, Mn{to) > 0, and 
Eq. (1191) then implies that all moments remain positive 
throughout the evolution, translating into positivity for 
/fc for all momenta at all times. 

Equation m closely resembles the Bobylev-Krook- 
Wu (BKW) equation derived almost four decades ago [12, 
[ 1 ^ in a famous study about homogeneous and isotropic 
solutions of the non-relativistic Boltzmann equation (for 


a review see 0 ) for Maxwell molecules. In fact, by 
dehning the time variable t = dt'/a^{t') to take into 

account the expansion [3^ , our equation (1191) for the mo¬ 
ments becomes identical to the BKW equations [l^ : 


drMmir) + Mmir) = 


TO-j-l 

3=0 


Mj{T)Mm-j{T). (20) 


This indicates that even though the underlying symme¬ 
tries of these physical systems are quite different (BKW’s 
are based on Galilean invariance with static conditions 
while ours are embedded in an expanding system), these 
systems are actually equivalent from a dynamical per¬ 
spective and evolve towards equilibrium in a universal 
manner. We note that our equations also reduce to 
BKW’s when a(i) = l, i.e. for a non-expanding metric 
(though special relativistic effects are still fully taken into 
account). 

4. Analytical solution of the moment equations. Given 
the close relation of our Eq. ([T^ with the BKW equation 
()20l) it is not surprising that it admits an exact analytic 
solution of the Krook and Wu type: 

Mmir) = IC{t)"^~^ [m - {m-l)IC{T)] (m > 0), (21) 

where/C(r) = 1 —exp(—r/6). The time evolution of the 
moments is shown in Fig. [T](a). One sees that low-order 
moments equilibrate more quickly than the higher-order 
ones that are needed to describe the high-momentum 
non-equilibrium tail of the distribution function. 

Given the exact form dSD of the moments Mm('r), 
the distribution function can be reconstructed as follows. 
One defines the function iF(T,u-k) = (u-k)^ 9{u-k) fk 
where u-k = k/a(T), expands its Fourier transform with 
respect to k in terms of the moments Mm{T), and then 
transforms it back. This yields the exact result jsij 

Mr) = Aexp(-^^^^) (22) 

- 4/C(r)-3 u-k / l-/C(r) Y 

_ /c4(t) r(r) V JCHt) )_ ■ 

To the best of our knowledge, this is the first analytic 
solution of the full Boltzmann equation for an expand¬ 
ing interacting gas. At the initial time io (corresponding 
to r = 0 and T(t)=Tq) one hnds the initial condition 
/fc(0) = ii (k/To) A exp[-4fc/(3To)] > 0. In Fig. [U (b) 
we plot the analytical solution for the ratio fk{T)/f^'^{T) 
as a function of u-k/T(T) = k/T^ for r = 0, 5, and 10. At 
r = 0 the high momentum tail is largely underpopulated 
whereas momentum modes in the range fc/Tg ^ 1.6—5 
are over-occupied relative to local equilibrium. As time 
evolves, high momentum modes are populated at the ex¬ 
pense of the over-populated moderate momentum region, 
in a process resembling an energy cascade. Note, how¬ 
ever, that when a{i) ^ one Hnds that limj^^,^ T(i) is 
finite and (1221) never assumes the equilibrium form. 
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k/To 


FIG. 1: (Color online) Time evolution of the moments 
(a) and of the ratio between the out-of-equilibrium solution 
122} and the equilibrium distribution (b). 


The moment equations (HU) can also be studied nu¬ 
merically for a variety of initial conditions for which no 
analytic solution is known. Here we briefly address the 
question whether the direct energy cascade seen in Fig.[T] 
is related to self-similar behavior characteristic of tur¬ 


bulence [ 33 . While a detailed study of self-similarity 
in this system is beyond the scope of this short paper, 
we note here that expressing the Boltzmann equation in 
terms of its moments may be quite useful in this context. 
For instance, one can show that non-thermal fixed points 
corresponding to exact self-similar solutions of the kind 


fkii) = a^{i) fs{a^{i) u-k), 


(23) 


where 7 and P are the scaling exponents and fs is the 
fixed point distribution (^ . are not compatible with 
Eq. 119}. In fact, using Eq. m one can see that the scal¬ 
ing exponents are necessarily 7 = 0 and /3 = 1, and that 
the Ansatz (123} for the distribution function (if assumed 
to be valid for all momenta) leads to moments 




1 


c 


s 

m 


a™+3(t) 


(24) 


(where c® = ^■m +2 Jg(^)) that have the same time 


dependence as the equilibrium moments. This implies 
that the corresponding normalized moments = 

Pm!Pm time independent. Since the conservation 
laws require that Mq = Mf = 1, one can then use HI} 
to show that t 1 for all m. This shows that there 

are no other true fixed points of the dynamics besides 
local thermal equilibrium. 

6. Conclusions. We derived from the relativistic Boltz¬ 
mann equation a nonlinear set of coupled moment equa¬ 
tions for a massless gas with constant cross section in a 
homogeneous, isotropically expanding spacetime. For a 
particular initial condition, we found that the moment 
equations can be solved exactly, thereby obtaining the 
first analytical solution of the Boltzmann equation with 
full non-linear collision term for an expanding system. 

The normalized moments of this expanding relativistic 
gas can be directly mapped onto the corresponding mo¬ 
ments for a homogeneous, static non-relativistic gas of 
Maxwell molecules. This happens even though the dis¬ 
tribution functions of these systems are not the same. 
This nontrivial correspondence suggests that exact self¬ 
similar non-equilibrium solutions may not exist even in 
rapidly expanding systems. 

The study performed here can be extended along sev¬ 
eral directions. One may consider different types of cross 
sections, or include nonzero particle masses to investi¬ 
gate bulk viscous effects. If possible, the generalization 
of the method presented here to exp anding systems with 
different symmetries [iS} lla, ISQ-lSSj that are relevant for 
the study of the quark-gluon plasma formed in heavy ion 
collisions would be particularly desirable. 
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